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Abstract. CO emission from the Sb(rs)I-II galaxy NGC 
7331 has been mapped in the J=2-l transition with a 21" 
beam over an area 3.5' by 1.3'. A relatively low contrast 
enhancement of molecular line emission occurs in a ring- 
like zone at a distance of approximately 3.5 kpc from the 
center; there is no evidence for a pronounced central hole. 
The ring is located at the edge of the region of rigid ro- 
tation and roughly coincides with an inhomogeneous ring 
of nonthermal radio continuum emission. It is well inside 
the radius of maximum rotational velocity. 

The intensities of the 492 GHz [CI] line and various 
^^CO and ^'^CO transitions observed towards the cen- 
ter and two outlying positions are modelled by multiple 
molecular gas components: low-density gas at a kinetic 
temperature Tkin ~ 10 K, and high-density gas at both 
2kin « 10 K and Tkin ~ 20 K. The molecular gas must be 
distributed in clumpy or filamentary form. The CO-to-H2 
conversion factor X applicable to the bulge is only half 
that applicable to the ring and beyond. The latter is still 
significantly lower than ^MiikyWay Molecular hydrogen is 
the dominant mass contributor to the interstellar medium 
in the bulge and in the ring. Far-infrared emission from 
dust peaks inside the ring at lOO/im (warm dust), and in 
the ring at 850/im (colder dust). Beyond the ring, neu- 
tral atomic hydrogen is dominant. Inferred total hydro- 
gen mass densities in the ring are about twice those in 
the bulge. Interstellar gas to dynamical mass ratios are of 
order 1% in the bulge, about 1.5% in the ring followed by 
a rise to 3%. The bulge gas may have originated in mass 
loss from bulge stars; in that case, the molecular ring is 
probably caused by a decrease in evacuation efficiency at 
the bulge outer edge. 

Key words: Galaxies - individual (NGC 7331) - ISM - 
spiral - structure; Radio lines - galaxies; ISM - molecules 



1. Introduction 

NGC 7331 is an isolated spiral galaxy of type Sb(rs)I-II 
with prominent dust lanes close to its centre (Kormendy & 

Send offprint requests to: F.P. Israel 



Table 1. NGC 7331 parameters 



Type-^ 


Sb(rs)FII 


Optical Centre: 




R.A. (1950)'' 
Decl.(1950)'' 


22''34'"47.7° 


34°09'35" 


Radio Centre : 




R.A. (1950)" 




Decl.(1950)'^ 


34°09'21" 


^LSR 


831 kms"^ 


Distance D*^ 


14.3 Mpc 


Inclination i-^ 


74.8° 


Position angle 


167° 


Luminosity Lg 


5.0 X 10^° Lb0 


Scale 


14.5 "/kpc 



Notes to Table 1: 

" RSA (Sandage & Tammann 1987) 

Dressel & Condon (1976) 
" Begeman (1987); Cowan et al. (1994) 

Corresponding to VhcI = 820 kms~^ (Begeman 1987) 
" TuUy (1988); corresponds to Ho = 75 kms"^ 

Hl-derived parameters from Begeman (1987) 
" Begeman (1987) rescaled to D = 14.3 Mpc 



Norman 1979; Sandage & Tammann 1987). Table 1 sum- 
marizes the relevant parameters of NGC 7331. A ringlike 
distribution of dust surrounding the bulge was suggested 
by Telesco et al. (1982). Such a distribution is also ap- 
parent in radio continuum maps (Cowan et al. 1994), but 
only vaguely in HI maps (Bosma 1978; Begeman 1987). 
A strong ring signature in CO emission was claimed by 
Young & ScoviUe (1982). 

Several galaxies are thought to host a molecular ring 
structure (see e.g. Braine et al. 1993), such as the ring in 
our own Galaxy discovered by ScoviUe & Solomon (1975) 
and the one in M 31 (Stark 1979; Dame et al. 1994; Koper 
1993). The latter serves to illustrate a dilemma commonly 
facing the interpretation of CO maps, especially those of 
highly inclined galaxies where rings are most easily dis- 
cerned: the conspicuous molecular structure may in fact 
consist of spiral arm segments that only in projection sug- 
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gest a ringlike structure. Two-dimensional CO mapping 
of NGC 7331 by von Linden et al. 1996) and Tosaki & 
Shioya (1997) support the latter interpretation. Although 
molecular rings have been identified in or claimed for other 
galaxies, the case of NGC 7331 is of interest because it re- 
sembles M 31 in being a large spiral galaxy of relatively 
early type, containing a prominent stellar bulge. It has also 
been claimed to have, as M 31, very little CO emission in- 
side its molecular ring (Young & Scoville 1982; Tosaki & 
Shioya, 1997). NGC 7331 even resembles M 31 in its high 
inclination (75° and 77° respectively - Arp & Kormendy 
1972, Sandage & Tammann 1987). Its radio structure is 
a stronger version of that of M 31 (Cowan et al. 1994). 
NGC 7331 also contains a clear, but patchy radio contin- 
uum ring. Inside the ring, little or no radio emission is 
found, except for a compact nuclear source. The luminos- 
ity of this source is 3-4 times that of Sgr A, and a thousand 
times stronger than the nucleus of M 31. It is associated 
with a nuclear X-ray source (Stockdale et al. 1998). Ring- 
like distributions of interstellar dust are also seen at mid- 
infrared (Smith 1998) and submillimeter (Bianchi et al. 
1998) wavelengths, but they are not nearly as evident at 
the far-infrared wavelengths inbetween (Smith & Harvey 
1996; Alton et al. 1998). 

An unusual characteristic of NGC 7331 is the rather 
low J=2-l/J=l-0 CO transitional ratio of 0.5-0.7 re- 
ported by Braine et al. (1993) and von Linden et al. (1996). 
This is quite different from most other galaxies observed 
in CO, where the two transitions are usually of similar 
strength. However, such low transitional ratios have also 
been found for individual dark clouds in the central parts 
of M 31 (Allen & Lequeux 1993). These were interpreted 
as evidence for very cold (Tkin < 5 K) and tenuous (n « 
100 cm~^) molecular clouds by Loinard et al. (1995), but 
Israel et al. (1998) showed that they are more likely caused 
by filamentary gas at temperatures Tkin « 10 K present 
at both low and high densities. Given the similarities be- 
tween NGC 7331 and M 31, it is of interest to investigate 
whether such a state of affairs also applies to the central 
region of NGC 7331. 

In this paper, we present a fully sampled map of NGC 
7331 in the J=2-l ^^CO transition over an area of 1.3' by 
3.5'. In addition, we have measured the first three ^^CO 
and ^■^CO transitions as well as the 492 GHz CI transi- 
tion towards the central region of the galaxy, allowing us 
to narrow down the permitted range of the apparently 
unusual physical conditions in the center. 



2. Observations 

Details relevant to the observations are listed in Table 2; 
the system temperatures given are the means for the re- 
spective runs. Observations in the J=l-0 transition were 
obtained with the IRAM 30 m telescope in service mode, 
at the optical and radio centre positions respectively, sepa- 



Table 2. Observations Log 



Transition 


Date 


Freq 


Tsys 


Beamsize 






(MM/YY) 


(GHz) 


(K) 


(") 
















J=l-0 


07/97 


115 


350 


21 


0.74 


J=2-l 


11/91 


230 


850 


21 


0.63 




Uo-lU/yz 




duu 


01 
Zi 


U.oo 




01/96 




500 


21 


0.69 




12/97 




260 


21 


0.69 


J=3-2 


05/93 


345 


720 


14 


0.53 




07/95 




900 


14 


0.58 














J=l-0 


07/97 


110 


200 


21 


0.74 


J=2-l 


09/93 


220 


660 


21 


0.63 




12/97 




320 


21 


0.69 


CI 












^Pi-^Po 


11/96 


492 


1900 


10 


0.53 



rated by Aa = 13.7", A5 = 14.0". The ^^CO observations 
were bracketed by the ^^CO observations. 

All other observations were carried out with the 15m 
James Clerk Maxwell Telescope (JCMT) on Mauna Kea 
(Hawaii) [|. Up to 1993, we used a 2048 channel AOS back- 
end covering a band of 500 MHz (650 km 8"^ at 230 GHz). 
After that year, the DAS digital autocorrelator system was 
used in bands of 500 and 750 MHz. Resulting spectra were 
binned to resolutions of 4 - 10 kms""'^, except for the 492 
GHz CI spectrum which was binned to 18 kms^^ on or- 
der to obtain a sufficiently high signal-to-noise ratio. Only 
linear baseline corrections were applied to the spectra. All 
spectra were scaled to a main-beam brightness tempera- 
ture, Tmb = T^/Vmh] relevant values for T^mb sue given in 
Table 2. Spectra are shown in Fig. 1 and summarized in 
Table 3. Position A (indicated in Figure 2) is that of the 
radio nucleus, the actual galaxy center. Along the minor 
axis, the beam just fills the space between the ring com- 
ponents. Positions B and C are close to the minor axis at 
deprojected radii R = 4.5 kpc and R = 8.2 kpc respec- 
tively; position B is centered on the CO ring. Because of 
the pronounced tilt of the galaxy, 21" circular observing 
beams sample elongated ellipses in the plane of the galaxy, 
covering a range of 5.55 kpc in the minor axis direction. 
The beams covering positions B and C overlap (see Figure 
2); positions A and B are essentially independent. 

For the J=2-l mapping observations, the integration 
time was typically 400 seconds per spectrum, on and off 

^ The James Clerk Maxwell Telescope is operated on a joint 
basis between the United Kingdom Particle Physics and Astro- 
physics Council (PPARC), the Netherlands Organisation for 
Scientific Research (NWO) and the National Research Council 
of Canada (NRC). 
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Fig. 1. Spectra observed towards positions A, B and C in 
NGC 7331 (see Sect. 2). Top row: J=l-0 CO; second row: J=2-l 

CO; third row: J=3 2 CO; bottom row: [CI]. Horizontal scale is 
LSR velocity in kms^^, vertical scale is in K. To convert to 
Tmb, multiply by 1.35 (J=l-0 CO), 1.45 (J=2-l CO), 1.72 
(J=3 2 CO) and 1.89 ([CI]), respectively. 

Note: CO is marked by solid lines, ^^CO multiplied by three is 
marked by dashed lines. 
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Table 3. CO and CI line intensities in NGC 7331 



Transition 


Resolution 


7mb 


/ Tmbdy 






( ) 


(mK) 


(K kms-i) 




A : Q = 


22:34:46.6; o 


= 34:09:21 


J =1-0 




21 


162 


26±2 




"CO 


21 


21 


3.9±0.4 


J=2-l 


i^CO 


21 


100 


14±2 




"CO 


21 


13 


2.5±0.6 


J=3-2 


i^CO 


14 


43 


7.5±1.3 






21 


49 


8.8±1.5 


^Pi-^Po 


CI 


10 


30 


1.9±0.3 




B : n = 


22:;i 1:17.7: d 


= :!!:()!): 


i") 


J =1-0 


12 CO 


■15" 


47 


7±1 






33'= 


170 


27±4 






01 


245 


25±2 






01 
Zi 


27 


3.3±0.4 


J=2—l 


CO 


Zi 


109 


12±2 






21 


13 


1.8±0.3 


J=3-2 


i^CO 




100 


10±2 






2r 




10±4 




C : a = 


22:34:48.9; <5 


= 34:09:34 


J =1-0 


i^CO 


21'* 




10±2 


J=2-l 


i^CO 


21 


36 


3.5±1 




"CO 


21 


5: 


0.6±0.3 


J=3-2 


i^CO 


14 


53 


2.5±0.5 






21"= 




3.5±1.5 



Notes to Table 3: 

a. Prom Young et al. (1995); b. Prom Elfhag et al. (1996); 
c. Extrapolated value; d. Prom von Linden et al. (1996) 



the source. After binning to a velocity resolution of 5 

kms""'^, the resulting r.m.s. noise and baseline deviations 
were of the order of 20 mK over most of the band. 

3. Results 

3.1. CO distribution 

In the J=l-0 and J=2-l transitions, the integrated inten- 
sities of both CO isotopes arc rather similar for positions 
A and B, indicating a fairly smooth central distribution 
of relatively strong CO. This is consistent with the large 
(50%) fraction of flux found to be missing by Tosaki & Sh- 
ioya (1997) in their J=l-0 CO interferometer map. The 
results presented here and by von Linden et al. (1996) 
arc, however, inconsistent with the folded major axis pro- 
file obtained by Young & Scoville (1982) and shown in 
more detail by Young et al. (1995), which exhibits both a 
pronounced lack of CO at the centre and a strong peak at 
a radial distance of 45" (3 kpc). From the data in Table 




+40 +20 -20 -40 +40 +20 -20 -40 
R.A. offset (arcsec) R.A. offset (arcsec) 



Fig. 2. Left: Digitized sky survey image of NGC 7331 with 
observed positions A, B and C (see Table 3) marked. Size 
of circle corresponds to 230 GHz beamsize. Right: Distri- 
bution of J=2-l *^C0 emission in NGC 7331, integrated 
over a velocity range Vlsr = 530 - 1130 km s~*. Con- 
tours are in steps of / TmhdV = 4 K km s~^. Observed 
positions A, B and C are indicated by circles. Table 1. 



2 and the results obtained by von Linden et al. (1996), it 

appears that the central integrated value given by Young 
et al. (1995) is too low by more than a factor of two. Thus, 
there is no significant central 'hole' in the distribution of 
CO emission, at least not on the scale of our 20" beam. 

A full-resolution contour map of the J=2-l CO in- 
tensity, integrated over the velocity range of 530 to 1130 
kms"-'^ is shown in Fig. 2. The map shows a good overall 
resemblance to the J=l-0 CO map obtained at slightly 
lower resolution by von Linden et al. (1996). This is also 
true for the major axis position-velocity diagram (not 
shown here). In Fig. 2, the elliptical outline of a low- 
contrast ring around the center can be discerned; peaks of 
CO emission occm at positions 30" north and 50" south. 
The latter two more or less correspond to the radial dis- 
tance of the molecular ring proposed by Young & Scoville 
(1982). The map covers the brightest part of the opti- 
cal image also whown in Figure 2 (for better images, see 
panel 40 in the atlas by Sandage & Bedke 1988). In this 
image, a large, overexposed bulge is surrounded by dust 
lanes and irregular spiral arms traced by Hll regions. The 
CO maxima at A6 = -50" and -|-30" fall on either side 
of the bulge, and the ring traces dusty spiral arms close 
to the bulge, especially on the western side. Most of the 
reddening of NGC 7331 occurs in this western spiral arm 
(Telesco ct al. 1982; Bianchi et al. 1998). On the eastern 
side of the CO map, faint emission due to a more distant 



F.P. Israel and F. Baas: Molecular gas in the bulge and ring of NGC 7331 



5 



major spiral arm is seen as well. The HI map obtained 
by Begeman (1987) at a very similar resolution shows an 
incomplete 'ring' of neutral hydrogen. The CO emission 
from the outlying spiral arm coincides with a relatively 
bright part of this HI 'ring'. Most of the CO emission is, 
however, well inside it and coincides with the radio con- 
tinuum ring mapped by Cowan et al. (1994). The main 
CO peaks are at the northern and southern extremities of 
the radio continuum ring. 

In the velocity-integrated single-dish CO map (Fig. 2), 
the ring is only weakly visible. Its presence is more clearly 
revealed in the interferometer map of Tosaki & Shioya 
(1997) and in Fig. 3. This figure shows the distribution 
of J=2-l ^^CO over the same region, but now integrated 
over velocity bins of 40kms~^ only. Between velocities 
Vlsr = 600 and 1000 kms-\ the maps show a double 
structure. The double-peak structure seen in Fig. 3 ex- 
tends over most of the part of NGC 7331 characterized 
by rigid rotation (cf. von Linden et al. 1996). A similar 
pattern, with the limitations imposed by interferometric 
techniques, is also evident in the channel maps published 
by Tosaki & Shioya (1997). 

3.2. Line ratios 

The ^^CO/^^CO isotopical ratios of 6 - 7 arc somewhat 
low compared to typical values around 10 found in most 
other external galaxies, but given the low IRAS /60/./100 
ratio of 0.3 (Rice et al. 1988) this is in line with the results 
obtained by Aalto et al. (1991). We find a ^^CO J=2- 
l/J=l-0 ratio of 0.54±0.10 in the center, consistent with 
the results obtained by Braine et al. (1993) and von Lin- 
den et al. (1996). Both J=l-0/J=2-l intensities and ratio 
decrease away from the center. The CO transitional ratios 
observed in NGC 7331 are rather different from those of 
most other galaxies, where the lower J transitions usually 
have similar velocity-integrated intensities (cf. Israel & 
van der Werf 1996). In contrast, CO intensities in the cen- 
tre of NGC 7331 decrease rapidly with increasing J level 
(Table 3). The J=2-l/J=l-0 ratios suggest subthermal 
excitation either at relatively low excitation temperatures 
or at very low column densities. However, the observed 
J=3-2/ J=2-l ratios of 0.6 or higher indicate the presence 
of a certain amount of warm molecular gas. Very low tem- 
peratures are also milikcly because they imply CO/ ^"^CO 
ratios substantially closer to unity than is observed. We 
also note, in Table 4, a similarity between the ratios ap- 
plicable to the D 478 cloud in the central parts of M 31 
and to the emission from NGC 7331, notwithstanding the 
factor of 400 difference in beam surface area. 

4. Analysis and discussion 

4.1. Modelling of observed line intensities 

We have modelled the observed intensities and their ra- 
tios by assuming the presence of two molecular gas com- 



ponents of different temperature and density, a relatively 
cold component dominating the J=l-0 emission and a 
warmer component becoming progressively more impor- 
tant in the higher transitions. We have used the radia- 
tive transfer models from the Leiden astrochemistry group 
(Jansen 1995; Jansen et al. 1994); we included a back- 
groimd radiation field of Tbg = 2.73 K. In these models 
kinetic temperature, molecular hydrogen density and CO 
respectively C column densities function as input param- 
eters. A further constraint is provided by the chemical 
models discussed by Van Dishoeck & Black (1988) which 
show a strong dependence of the N{C)/N{CO) column 
density ratio around molecular hydrogen column densities 
of about 10^1 cm-2. Above N{ H2) = 2 x lO^i, virtually aU 
carbon is in CO, whereas below N{Yi2) = 2 x 10^°cm~^ 
virtually all carbon is in C. The CI and CO intensities ob- 
served in position A, together with the sensitivity of the 
C/CO ratio to H2 column density provide rather strin- 
gent constraints on the models acceptable for at least the 
center of NGC 7331. 

Although the available data do not allow a precise 
and unique determination of the physical condition of 
the molecular gas in NGC 7331, they serve well to con- 
strain the parameter space of possible solutions. For in- 
stance, cold component CO column densities cannot ex- 
ceed N{CO) « 10^^ cm^^ and warm component kinetic 
temperatures cannot be lower than 20 K without con- 
fiicting with observed -'^^CO and -'^^CO ratios. In Table 5 
we list a few representative models. Note that the CO 
column densities listed in Table 5 are those of a sin- 
gle model cloud. We assume a homogeneous population 
of such model clouds, non-shadowing in position- velocity 
space, so that the actual galaxy beam-averaged column 
density is the sum of the model cloud column densities in 
the beam multiplied by their beam filling factor. 

In models 4 through 7 we vary the cold component 
input parameters, and assume a warm component of 30 
K and density n{B.2) = 1000 cm~"^. In models 2 and 3, 
we have changed the warm component temperature to 40 
K and 20 K, and in model 1 we assume a more com- 
plex situation. In addition to the warm component, the 
cold component itself is structured into a high-density 
and a low-density contributor. It is a scaled version of 
the single-temperature, dual-density model (Tkin = 10 K; 
n(H2) = lOOcm-3 and 3000 cm'^) applied to the M 31- 
D 478 cloud complex by Israel et al. (1998). Although 
a warm component must be included, its nature is un- 
clear. Given the large linear beamsize (1.5 x 5.5 kpc) in 
the plane of the galaxy, this warm component may rep- 
resent discrete, starforming cloud complexes at some dis- 
tance from the center. In all models, J=l-0 CO intensi- 
ties are dominated by emission from the cold component 
with contributions of about 75%, 70% and 85% for posi- 
tions A, B and C respectively. In contrast, the J=3-2 CO 
intensities are all dominated by emission from the warm 
component. For the optically thin ^^CO transitions the 
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50 -50 

R.A. offset(arcsec) from 22 34 47.26 



Fig. 3. J~2-l ^^CO channel maps of NGC 7331. Emission is integrated over velocity bins of 40 kms ^; central 
velocities are marked in the panels. Contours are in steps of / T^nbdV = 2 K km s~^. 



Table 4. Integrated line ratios in the centre of NGC 7331 



Transitions 




pos. A 


pos B. 


pos. C 


D 478" 




(2-l)/(l-0) 


0.54±0.10 


0.51±0.10 


0.35±0.11 


0.42±0.08 




(3-2)7(2-1) 


0.63±0.15 


0.78±0.35 


1.0±0.41 


0.33±0.10 


"CO 


(2-1)7(1-0) 


0.64±0.2 


0.55±0.1 




0.45±0.15 


^^co/"co 


(1-0) 


6.7±1.1 


7.8±1.2 




8.8±1.8 




(2-1) 


5.6±1.6 


7.2±1.6 


5.8±3.3 


8.4±1.7 


CI/C0(2-1) 




0.21±0.07 






0.24±0.05 



Note: a. Dark cloud in M 31; see Allen et al. (1995); Loinard & Allen (1998); Israel et al. (1998). 



Table 5. Model parameters for NGC 7331 



Model 




Cold Component 


Warm Component 




Kinetic 


Gas 


CO Column 


Kinetic 


Gas 




Temperature 


Density 


Density 


Temperature 


Density 




Tkin 


n(H2) 


N(C0)7dV 
(10^^cm-2(Kkms-^)-^) 


Tkin 


n(H2) 




(K) 


(cnr") 


(K) 


(cm-3) 


1 


IO7IO" 


10073000 


0.777.0 


30 


3000 


2 


10 


300 


1.0 


40 


1000 


3 


10 


300 


1.0 


30 


1000 


4 


10 


300 


1.0 


20 


1000 


5 


10 


100 


1.0 


30 


1000 


6 


10 


300 


0.1 


30 


1000 


7 


10 


100 


0.1 


30 


1000 



Note: " Cold component assumed to consist of both low and high density gas as in D 478; see Israel et al. (1998). 
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Table 6. CO and C in NGC 7331 



Model 


Predicted Beam- Averaged 


Conversion 


Mass per Beam 


Face- 


■on 




CI Intensity Column Densities 


Factor X 






Mass Density 




fT^bdy N{CO) N{C) iV(H2) 




M(H2) 




a(H2) 






(Kkms-i) (102°cm-2) (10^° cm-VKkms-^) 


(10^ 


Mo) 


(Mo/pc-^) 




rOSlLlOll i\ H /JVC — -L / UU, iV \ll 1 ) — OX lU 


_2 

cm 










1 


2 0.011 0.006 12 


0.4 


3 


5 


5 


8 


2 


4 0.012 0.040 42 


1.6 


11 


15 


17 


24 


3 


4 0.009 0.041 40 


1.5 


10 


14 


16 


23 


4 


4 0.008 0.042 40 


1.5 


10 


14 


16 


23 


5 


4 0.009 0.051 48 


1.9 


13 


18 


20 


28 


6 


8 0.005 0.045 39 


1.5 


10 


15 


16 


23 


7 


9 0.008 0.089 80 


3.1 


21 


29 


33 


46 




Position B; Nh/Nc = 2300; N{HI) = 9 x 10^° 


cm^'^ 










1 


2.5 0.011 0.007 16 


0.6 


4 


7 


7 


12 


3 


4 0.006 0.030 37 


1.5 


10 


15 


15 


23 




Position C; Nh/Nc = 3000; N{HI) = 9 x 10^° 


cm~^ 










1 


0.7 0.004 0.003 8 


0.8 


2 


4 


3 


7 


3 


1.3 0.002 0.013 19 


1.9 


5 


8 


8 


13 



situation is loss cloarciit; if wo assume an intrinsic isotopic 
ratio of 100, emssion from the cold component contributes 
about 25% to the J=l-0 emission from positions A and B, 
whereas this fraction increases to about 45% if we assume 
an isotopic ratio of 50. 

4-2. CO and C column densities 

In order to relate neutral carbon and carbon monoxide 
column densities to that of molecular hydrogen, we have 
used [C]/[H] gas-phase abundance ratios estimated from 
the [0]/[H] abundance. From the data tabulated by Zarit- 
sky et al. (1994) we determined for the central beam (po- 
sition A) 12 + log (0/H) ^ 9.2, i.e. [0]/[H] = 1.5 x IQ-^. 
Although high, such an oxygen abundance is normal for 
galaxy centers (Garnett et al. 1997; van Zee et al. 1998). 
Using results given by Garnett et al. (1999), notably their 
Figures 4 and 6, we arrive at an estimated carbon abun- 
dance [C]/[H] = 2±1 X 10~^. As a significant fraction 
of all carbon will be tied up in dust particlc^s, and not 
be available in the gas-phase, we adopt a fractional cor- 
rection factor (5c = 0.33. Neglecting contributions by e.g. 
^■^CO and ionized carbon, we thus find A^h = [2A^(i?2) + 
N{HI)] K. 1700 [N{CO) + N{C)] with a factor of two 
uncertainty in the numerical factor. Similarly, we find for 
the off-center positions B and C numerical factors of 2300 
and 3000. The beam-averaged column densities in Table 
6 have been obtained by scaling the model cloud column 
densities by the ratio of actual observed CO intensity to 
predicted model CO intensity. 

The results of our model calculations are given in Ta- 
ble 6. In the table we give the predicted [CI] intensity 
Ici, which can be verified observationally, the calculated 



beam-averaged cohimn densities for both CO and C, the 
H2 column densities derived from these using the A''^ /Nc 
ratios and N{HI) values given, as well as the implied CO 
to H2 conversion factor X = N{Il2)/ Ico- The neutral 
carbon intensities I{CI) were calculated under the as- 
sumption that a significant fraction (0.6-0.7) of the total 
atomic carbon column density is ionized and present in 
the form of [CII]. Changes in the input CO column den- 
sities do not strongly affect the resultant C column den- 
sity: a substantially higher N{CO), for instance, implies a 
lower N{C)/N(CO) ratio, yielding a relatively unchanged 
N{C). We have also performed the calculations for a ratio 
of 50. Generally, the ratio of 100 provides a better fit to 
the ^"^CO data than the ratio of 50. As the end results 
for the two sets are moreover very similar, we have not 
included the latter in the table. The results of all models 
are given for position A, where we have also measured the 
[CI] intensity in a 10" beam. However, the models apply 
to measurements in the 21" beam observed or synthesized 
for CO. If atomic carbon is at a minimum in the center, 
[CI] intensities in a twice larger beam may be somewhat 
higher, perhaps by as much as 40%. Table 6 shows that 
model 1 yields a very good fit, whereas models 2 through 
5 are marginally possible and models 6 and 7 are ruled 
out. Model 1 is not unique; various other combinations 
of somewhat different kinetic temperatures for both cold 
and warm gas and somewhat different densities, yield very 
similar results. 

As models 6 and 7 are ruled out for position A and 
models 2 through 5 yield almost identical final results, we 
present only models 1 and 3 for positions B and C. The re- 
sults for position C have relatively large uncertainties due 
to the weakness of its emission, and the lack of a J=l-0 
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'^"^CO measurement. The results are not greatly different 
from those obtained at position A. Column densities de- 
crease, and X factors increase somewhat with radius. At 
Tkin = 10 K, the [CI] transition at 809 GHz has 

negUgible intensity, but this becomes comparable to the 
3Pi-^Po 492 GHz transition at Tki„ = 30 K. The presence 
of the warm component can therefore be verified by future 
observations of the 809 GHz [CI] transition, for which we 
predict an intensity of 15-30% of the 492 GHz intensity. 
For the [CII] emission we expect intensities of the order of 
5 X 10~^ erg s""'^ cm~^ sr"-"-. 

Beam-averaged neutral carbon to carbon monoxide 
column density ratios are N{C)/N{CO) = 0.65±0.1 and 
N{C)/N{CO) = 5.5±1.0 for models 1 and 3 respectively. 
The former is close to the typical values 0.2-0.5 found for 
M 82, NGC 253 and M 83 (White ct al. 1994; Israel ct al. 
1995; Stutzki et al. 1997; Petitpas & Wilson 1998), but the 
latter is much higher and is only matched by the corre- 
sponding ratio of 3-6 found in Galactic translucent clouds 
(Stark & van Dishoeck 1994). 

^.3. Molecular hydrogen and the I{CO) to N{'R2) ratio 

Although any explanation of the observed CO intensities 

requires the presence of both cold and smaller amounts 
of warm molecular gas in NGC 7331, the range of admis- 
sible parameters is not fully constrained. The [CI] inten- 
sity observed towards the center of the galaxy, however, 
strongly suggests a complex physical environment of the 
sort represented by model 1. This model is characterized 
by cold molecular gas (typical temperature Tki„ = 10 K) 
present at both high and low volume densities (typically 
of order a few hundred and a few thousand per cc respec- 
tively), in addition to a warmer component (temperature 
Ikin > 20-ftr) of high density. This is probably a simplifi- 
cation: in reality a range of densities and temperatures is 
likely to be present. As the large linear beamsizc (1.5 kpc 
along the major axis, 5.5 kpc along the minor axis) only 
provides results averaged over a large radial range, the 
spatial distributions of the cold and the warm gas within 
the beam may well be different. Both kinetic temperature 
and mean molecular gas density in the centre of NGC 7331 
are typically an order of magnitude below the values found 
in later-type starburst galaxies such as NGC 253 and M 82 
(Israel et al. 1995; Wall et al. 1991). 

Our models suggest that a large fraction of the CO 
emission originates from cold gas of low column density. 
A smaller fraction originates in much denser gas, partly 
at higher temperatures. The cold gas in the center of 
NGC 7331 appears to be similar to that in cloud com- 
plexes such as D 478 in M 31 dark and the Taurus- Auriga 
complex in the Milky Way; most likely, it is highly frag- 
mented and filamentary (Israel et al. 1998). The warm gas 
may be heated by the nucleus and by luminous stars in 
the inner spiral arms and the 'molecular ring'. The pres- 
ence of energetic photons in the inner part of NGC 7331 



is betrayed by Hq:+[NII] emission (see Fig. 5 by Smith & 
Harvey 1996) and more directly by significant UV emis- 
sion (Wesselius et al. 1982) unlikely to be dominated by 
the spiral arms because of their high dust content (Bianchi 
et al. 1998). 

The evaluation of the models in Tables 5 and 6 as- 
sumes a radiation field /uv ~ 1, corresponding to /looo = 
4.5 X 10~* photons s~^ cm~^ which is consistent with 
the longer- wavelength UV data by Wesselius et al. (1982). 
Such a low radiation field density is also indicated by the 
strength of the 7.7 and 11.3 /xm dust emission features 
(Smith 1998). The beam-averaged column densities in Ta- 
ble 6 are relatively insensitive to changes in the assumed 
/uV) because in the cold diffuse gas, most carbon is already 
in C rather than in CO, whereas the much smaller filling 
factor of the dense gas greatly reduces the effect of changes 
in the N(C)/N(CO) ratio on the beam-averaged neutral 
carbon column density. Moreover, we expect only limited 
variation (by a factor of 2-3) in the radiation field density 
over at least the inner 5 kpc because of the smooth distri- 
bution of Ha emission as well as the far-infrared emission 
between 50/im and 200/im (Smith & Harvey 1996; Alton 
et al. 1998). The quiescence of the spiral arms is illustrated 
by the strong excess of 450 and 850/im emission from cold 
dust (Bianchi et al. 1998). Thus, the spiral arms contain a 
relatively large amount of cold dust especially in compar- 
ison with the central region. We are therefore confident of 
the derived iV(H2) values in Table 6. 

As C and O abundances in NGC 7331 are 2-5 times 
higher than those in the Solar Neighbourhood, and radi- 
ation fields are not particularly intense, we expect CO in 
NGC 7331 to be relatively well-shielded, so that the CO 
to H2 conversion factor X should be lower than that in 
the Solar Neighbourhood, i.e. fewer H2 molecules per unit 
CO intensity. Indeed we find values of X lower than the 
value of 2 X 10^°cm~^/Kkms~^ assumed for the Milky 
Way, which can also be compared to the relationship be- 
tween X, radiation field intensity and metallicity found by 
Israel (1997). For positions A, B and C we take [0]/[H] 
abundance ratios of 1.5, 1.3 and 1.1 in units of 10~^ re- 
spectively (Zaritzky et al. 1994). From high-resolution far- 
infrared surface brightnesses (Smith & Harvey 1996), HI 
column densities (Begeman (1987) and eqn. (3b) from Is- 
rael (1997), we predict values X = 0.8, 0.7 and 1.0 in 
units of 10^'^ cm~^/(Kkms~^)~^ for positions A, B and 
C respectively. These are very close to the results from 
the preferred model 1, and a factor of two or more below 
the results for the other models. Neglect of the radiation 
field term in Israel's (1997) eqn. (3b), i.e. use of his eqn. 
(4) predicts in the same units X = 0.15, 0.25 and 0.4 for 
positions A, B and C, i.e. much lower than any of the 
model results. We conclude that the low values of X in 
the preferred model 1 are in good agreement with both 
the high abundances in NGC 7331 and the relationship 
between X , radiation field intensity and metallicity found 
by Israel (1997). 
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With respect to the value of X derived for position 
A it should be noted that the large linear beamsize in- 
cludes both the center of NGC 7331 and more outlying 
regions along the minor axis. If the latter were to be char- 
acterized by an X value closer to that of position B, the 
actual central X value would be significantly lower. For 
instance, if we assign X = 0.6 x 10^° cm~^/ K kms""'^ to 
the outer half of the CO emission, the inner half would 
have X = 0.2 x 10^'^cm~^/Kkms~^, an order of magni- 
tude less than the Milky Way value, and well below what 
is suggested by the high metallicity. Such a very low value 
would, however, not be unexpected. For the Milky Way 
centre, Sodroski et al. (1995) conclude to an X-factor 3 - 
10 times smaller than the 'standard' Galactic value. The 
COBE Galactic Centre data presented by Bennett et al. 
(1994) imply lower CO transition ratios somewhat similar 
to those in NGC 7331. 

Another way of verifying the derived H2 column den- 
sities is provided by the submmillimeter observations pre- 
sented by Bianchi et al. (1998). For Fgso^m = 50 mJy and 
T = 20±3 K in a 30" x 40" beam (Bianchi et al. 1998), we 
derive a beam-averaged = 6.6 (-1.2, +2.2). Further- 
more assuming that the dust to gas ratio is proportional 
to metallicity, we modify the Galactic relation between to- 
tal hydrogen column density and visual extinction (Bohlin 
et al. 1978) to Nn = 0.6 x lO^^Ay cm-^. This imphes a 
column density A''(H2) = 2(-0.4, -1-0.6) x 10^^ cm'^ (cor- 
responding to X = 0.75(-0.15,+0.25) x 10^"). The simi- 
larly obtained result for position B is slightly lower. These 
results are thus in rather good agreement, given the var- 
ious uncertainties, with 7V(H2) = 1.2 — 1.6 x 10^^ cm^^ 
and X = OA - 0.6 x 10^" found for positions A and B 
using model 1. 

Comparison of the models and the observations allows 
us to draw some general conclusions on the distribution 
of molecular hydrogen in NGC 7331. In model 1, rela- 
tive amounts of cold/tenuous, cold/dense and warm/dense 
molecular hydrogen gas are 45%, 30% and 25% for posi- 
tions A and B. The results for position C seem to indicate 
a somewhat higher contribution by warm molecular gas. 
Using the beam-averaged H2 column densities and the 
model H2 volume densities, we find that the average line 
of sight within the beam contains cold/tenuous II2 over 
about 2 pc (Model 1, positions A and B) to 0.7 pc (Model 
1, position C). Both the cold and the warm dense com- 
ponent have average line-of-sight extents a factor of 50 
lower. However, the observed CO temperatures are much 
lower than the model excitation temperatures, indicating 
small beam-filling factors for the molecular material. As- 
suming individual lines of sight within the beam to be 
either empty, or homogeneously filled with molecular gas, 
we find for those line of sights that do contain molecu- 
lar gas extents of about 20 pc (cold tenons gas), 2.5 pc 
(cold dense gas) and 25 pc (warm dense gas); these num- 
bers are indicative of the max;imum source size that can 
be expected. 



Application of model 3 yields somewhat different 
results. Here, most of the molecular gas is in the 
cold/tenuous form rather than in the warm/dense phase: 
only 17% warm gas is required at position A, and about 
5% at positions B and C. Average line of sight extents are 
3.5 pc for cold gas at positions A and B, and half that at 
position C. After correction for beam filling, we find lines 
of sight extents of typically 115 pc for the cold molecular 
gas and 15% or less of that for warm gas. Only at position 
C a more uncertain extent of 20 pc is obtained. 

The derived line of sight extents are much smaller than 
the length of the line of sight traversing the galaxy, which 
is about four times its thickness. Although the latter is 
not known, this length can be estimated at well in excess 
of a kiloparsec. Thus, only a small fraction of the volume 
sampled by the beam at each of the analyzed positions 
is filled with molecular material. This material is highly 
clumped or distributed in filamentary form. 

4-4 ■ Radial distribution of molecular gas 

In the case of highly-inclined ring structures, major-axis 
position-velocity diagrams may give a misleading impres- 
sion of the actual radial distribution of emitting material, 
becaiise at the tangential points substantially longer lines 
of sight contribute to the emission. To determine the ac- 
tual distribution of CO as a function of radial distance 
from the centre spectra, we have fitted an inclined axisym- 
metric disk model to the data in the velocity-integrated 
map (Fig. 3) by applying the Richardson-Lucy iterative 
scheme (Lucy 1974). In principle, with a priori knowl- 
edge of the (CO) velocity field, this technique can also be 
used to obtain radial distributions with a spatial resolution 
higher than that of the observing beam (Scoville, Young 
& Lucy 1983). Fig. 4 shows the fitted radial distribution 
of the velocity-integrated J=2-l CO emission. 

The fitted profile corresponds to the face- on, radial dis- 
tribution of /'(CO)o = / {T*{CO)dV)o. The CO luminos- 
ity starts at I'{CO)o = 4 K kms~^ in the center, reaches 
a minimum at i? « 1.75 kpc, and the reaches a maximum 
at R « 3.5 kpc after which it drops smoothly to Iqq = 
2 K kms~^. The ring-to-disk intensity contrast ratio is 
about 0.6. The molecular 'ring' is clearly discernible, but 
it does not dominate the CO distribution in the galaxy. 
Both the major-axis CO distribution and the fitted radial 
CO profile are different from those of M 31, where most 
of the CO is found farther out in the spiral arm 'ring' at 
i?=9 kpc, and very little CO occurs at the centre (Dame 
et al. 1994). 

In order to determine the radial distribution of inter- 
stellar gas in NGC 7331, we have converted the CO radial 
profile to a radial distribution of H2 mass densities, using 
the J=l-0/J=2-l CO ratios and X values from Tables 
4 and 6, and combined these with the radial HI profile 
published by Begeman (1987). Both the radial H2 and 
HI profiles are also shown in Figure 4. In the inner 4.5 



10 



F.P. Israel and F. Baas: Molecular gas in the bulge and ring of NGC 7331 



kpc, the H2 mass dominates that of HI by about 40%. 
Beyond this, HI becomes increasingly dominant. The ra- 
dial distribution of molecular gas reaches its peak (at R 
=3.1 kpc) well before that of HI (at i? « 10 kpc). Al- 
though the radial J=2-l CO profile exhibits a relatively 
low contrast between the ring feature and the underlying 
disk emission, the radially increasing transitional ratios 
and X-values serve to enhance the contrast in H2. The 
central region is not empty, but the H2 mass density in- 
side the ring is only 75% of that in the ring; due to the 
lack of HI in the center, the relative mass density of all 
hydrogen is even lower with 55% of the ring value. The 
facc-on H2 mass distribution increases from cru^ = 6 Mq 
pc^^ to 8 M0 pc~^ at i? = 3.1 kpc. The total hydrogen 
radial mass-density distribution increases from a central 
value (THI+H2 = 8 Mq pc~^ to (Thi+H2 = 14 Mq pc~^ at R 
= 3.5 kpc and then drops slowly. The gaseous fraction (in- 
cluding helium) of the total mass was estimated from the 
rotation curves given by Rubin et al. (1965) and Begeman 
(1987), assuming a spherical bulge and circular velocities. 
Inside the ring, the gas-to-total mass ratio Mgas/Mdyn is 
about 1%. In the ring, it rises to 1.5%, and then slowly 
climbs to 3%. From Begeman's (1987) data, neglecting H2, 
we find in comparison a global ratio Mgas/Mdyn of 3.2%. 
Even with the dominant contribution by H2, the gas in 
the inner part of NGC 7331 is only a minute fraction of 
the total mass. 

Finally, it is of interest to compare the radial gas dis- 
tributions to the radial far-infrared profiles representing 
interstellar dust. The 100/um far-infrared profile (Smith & 
Harvey 1996) is rather flat in the center, reaching a very 
minor maximum at about 2 kpc, after which a steep de- 
cline sets in. The 450/Lt and 850/um profiles show radially 
increasing intensities that reach a peak at about 3 kpc, the 
850/im profile peaking a little farther out than the 450 /xm 
profile (Bianchi et al. 1998). The H2 mass-density peaks at 
R = 3.1 kpc, just about coincident with the 85G/im maxi- 
mum. At the location of the molecular 'ring' peak, lOO/xm 
intensities are about 80% of those in the center. However, 
the total gas mass-density peaks slightly farther out at R 
= 3.8 kpc. Thus, emission from warm dust peaks well in- 
side the molecular ring, and both molecular gas and cold 
dust peak inside the radius of highest gas mass-density. 
Beyond the ring peak, the decline of the 450/x and 850/xm 
emission more or less follows that of the H2 mass density 
profile. We conclude that the interstellar dust is hottest 
in the central region, where gas mass-densities are lowest. 
The mean dust temperature defined by the far-infrared 
ratios smoothly decreases from the center reaching a shal- 
low minimum at about i? = 3 kpc, i.e. at the H2 peak, 
beyond which it appears to increase slightly. 

4-5. Origin of bulge molecular gas and the ring 

Various mechanisms for the occurence of ring morpholo- 
gies have been suggested in the literature. Interaction be- 
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Fig. 4. Deprojected radial profiles. Bottom: Face-on radial 
distribution of CO emission obtained with the Richardson- 
Lucy scheme (see section 4.4). Vertical axis is / TmbdV as 
would be observed perpendicular to the galaxy plane. Top: 
Face-on mass-densities of H2 and HI, in units of M© pc~^. 
HI data were taken from Begeman (1987). 

tween magnetic fields and the gas distribution were pro- 
posed and tested by Battaner et al. (1988). Inner Lind- 
blad resonances may be responsible (Kormendy & Norman 
1979), and a ringlike feature may also result from evacu- 
ation of gas from the central regions by stellar winds (cf. 
Faber & Gallagher 1976; Soifer et al. 1986, Muhoz-Tuhon 
& Beckman 1988). Finally, as Young & Scoville (1982) sug- 
gested for NGC 7331, a ringlikc appearance may also be 
caused by the nuclear bulge having used up the originally 
present molecular gas in the center. 

In NGC 7331, the solid-body rotation curve rises 
rapidly out to R = 3.5 kpc after which it flattens and 
reaches a broad maximum at J? = 6 kpc (see also von Lin- 
den et al. 1996). The ring is thus located just at the radius 
where rigid rotation is lost, but well within the radius of 
maximum rotational velocity (iiring = 0.6 i?vmax)- This is 
unlike M 31, where the molecular ring is found at a radius 
twice that of peak rotational velocity (Brinks & Burton 
1984; Dame et al. 1993). Both the molecular ring and the 
boundary of the solid-body rotation region are also just 
at the radius at which the light of the disk becomes domi- 
nant over that of the bulge (cf. Begeman 1987). Although 
the radius of the mostly nonthermal radio continuum ring 
radius is more difficult to determine because of its inhomo- 
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geneous structure, it appears to coincide more or less with 
the molecular ring, also well inside the radius of maximum 
rotational velocity (Cowan et al. 1994). 

On the same reasoning as used by Yoimg & Scoville 
(1982), we may rule out the presence of an inner Lind- 
blad resonance in NGC 7331 as an explanation for the 
observed molecular ring structure, because an ILR can 
only occur well outside the region of solid-body rotation 
(e.g. Kormendy & Norman 1979). Our observations clearly 
show that there is no pronounced CO hole in the centre 
of NGC 7331, although the derived distributions of both 
B.2 and total interstellar gas do show a significant central 
depression, von Linden et al. (1996) have suggested that 
the ringlike distribution in NGC 7331 is caused by the dy- 
namical action of a weak central bar. However, their simu- 
lations yield both a ring more massive than observed, and 
center more devoid of gas than observed, casting doubts 
on the proposed central bar. 

In the case of M 31, Soifer et al. (1986) suggest that 
the amount of interstellar matter observed in the center of 
M 31 could have accumulated from late- type stellar mass 
loss in the bulge, and is kept low by continuous gas re- 
moval by supernova explosions and star formation. Could 
this also be the case in NGC 7331? It has been suggested 
by Prada et al. (1996) that the bulge of NGC 7331 is 
counter-rotating. Since the bulge gas is rotating in the 
normal sense, this would seem to preclude a stellar ori- 
gin for the gas. However, spectroscopy by Mediavilla et 
al. (1997) and Bottema (1999) does not confirm the sug- 
gested counter-rotation. The total mass of the interstellar 
gas inside R = 2 kpc is 1.6 x 10** Mq . According to the rea- 
soning outlined by Soifer et al. (1986), stellar mass loss in 
the bulge would accumulate this amount in 4 x 10^ years. 
Removal of the same amount of material from the bulge 
requires the energy output of 6/nx 10^ Type I supernovae, 
n being the fraction of energy available for the accelera- 
tion of interstellar material. For a Type I SN rate of 4 x 
10~^^ Lb yr"^ (Lang, 1992) the timescale for removal is 
6.5/n X 10^ years. Only if the fraction of supernova en- 
ergy actually available for removal exceeds 17%, will the 
interstellar gas be evacuated from the bulge faster than 
bulge stars can manufacture it. More generally, with the 
assumptions from Soifer et al. (1986), the ratio of evacu- 
ation to deposition timescales is tg/td, = 3.75 x lQ~^v^/n. 
The tabulation of NGC 7331 rotation velocities by Bege- 
man (1987) then suggests relatively efficient evacuation in 
the inner 1 kpc {te/td = 0.01/n — 0.05/n), and much less 
efficient evacuation at the edge of the bulge (i? = 4 kpc; 
te/td = 0.24/n). The radial decrease of the ratio of far- 
infrared emission to H2 mass-density found above may be 
related to this finding. It thus appears that the relatively 
small amounts of interstellar gas in the bulge of NGC 7331 
(Mgas/Mdyn ~ 0.01) also may well be the result of mass 
loss from the bulge stars themselves, rather than the re- 
sult from a net inflow of molecular material from greater 
radii. 



5. Conclusions 

1. Analysis of the J=2 1 ^^CO distribution and kinemat- 
ics shows the presence of enhanced molecular emission 
in a ringlike zone in NGC 7331, peaking at a radial 
distance of 3.5 kpc with a width of about 2 kpc. At 
R = 3.5 kpc, the velocity-integrated CO intensity of 
the ring itself is about 0.6 times that of the underlying 
more smoothly distributed CO emission that fills the 
entire bulge of NGC 7331. 

2. The velocity-integrated CO intensities in the center of 
NGC 7331 decrease strongly with increasing rotational 
level. The intensities in the J=l-0, J=2-l, J=3-2 
transitions are in the ratio of 1.0 : 0.55 : 0.35 respec- 
tively. The observed ^^CO/^'^CO isotopic ratios are 6.7 
and 5.6 in the J=l-0 and J=2-l transitions respec- 
tively. Positions at larger radial distances have similar 
ratios, albeit with somewhat stronger J=3-2 CO emis- 
sion, and weaker ^^CO emission. Weak [CI] emission 
was detected from the center. 

3. Modelling of the observed line ratios suggest a multi- 
component molecular medium. Gas with a kinetic tem- 
perature of about 10 K appears to be present at both 
low and high densities. At high densities, a warmer 
component with a kinetic temperature of 20 K or more 
is also present within the observing beams. The gas 
is probably distributed in a clumpy and filamentary 
form. 

4. Assuming a [C]/[H] abundance ratio of the order of 1-2 

xlO"'', the mean CO-to-H2 conversion factor is Xntssi 
= 4 X 10^^ cm~^ in the bulge region, and double that 
value in the ring and beyond. These values are well be- 
low those found in the Solar Neighbourhood, but they 
are consistent with the high metallicity of NGC 7331 
and with submillimeter dust observations. 

5. In the bulge, interstellar gas (HI -I- H2 -I- He) mass 
densities, projected onto the plane of the galaxy, are 
of the order of 11 M0 cm^^. In the ring itself, now 
properly placed at i? = 3.1 kpc, the gas mass density 
is almost twice as high. Within the ring, the interstel- 
lar gas mass is dominated by the molecular hydrogen 
contribution. Gas to total (dynamical) mass ratios are 
about 1 % in the center and about 1.5 % in the ring. 

6. The molecular ring coincides more or less with the 
mostly nonthermal radio continuum ring and the 850 
/im ring representing emission from cold dust. Emis- 
sion from warmer dust in the 100/xm wavelength range 
peaks well inside the molecular ring; dust temperatures 
appear to be decreasing with radius reaching a minin- 
mum in the ring. The radial distribution of HI reaches 
it maximum well beyond the molecular ring. 

7. The molecular ring is well inside the radius of peak 
rotational velocity. Its maximum is just at the edge 
of the region of solid-body rotation, and just at the 
radius where disk light becomes dominant over bulge 
light. The ring is not associated with an inner Lind- 
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blad resonance. The molecular gas inside the ring may 
have originated from mass loss by late type stars in 
the bulge. If this is the case, the ring is probably the 
result of wind-driven gas removal from the center. 
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